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Розглянуті похибки, які виникають при перетворенні типу «область - область». Запропоновано метод 
перетворення областей зображень, який базується на перетворенні скелетів областей. Розроблений 
алгоритм визначення похибки перетворення типу «область -- область». Розроблені алгоритми програмно 
реалізовані в середовищі Вогіапа Дербі 7.0. Досліджені похибки перетворення типу «область -- область» в 
автоматизованих системах аналізу цитологічних зображень. 
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Рассмотрень погрешности, возникающие при преобразований типа «область - область». Предложен 
метод преобразования областей изображений, основанньїй на преобразований скелетов областей. Разработан 
алгоритм определения погрешности преобразования типа «область -- область». Разработаннье алгоритму, 
программно реализованьєе в среде Вогіапа Де!рбі 7.0. Исследовань погрешности преобразования типа 
«область - область» в автоматизированньх системах анализа цитологических изображений. 
Ключевье слова: погрешности преобразования типа «область -- область», скелеть |, 

цитологические изображения, автоматизированньге системь. 


Вступ 


Сучасні автоматизовані засоби опрацювання біомедичних зображень (АЗОБМЗ) 
складаються в основному з таких модулів: модуля вводу зображень, який представ- 
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лений світловим мікроскопом 1 відеокамерою, персональним комп'ютером, котрий 
забезпечує опрацювання зображень і принтером для документування результатів 
аналізу на папері (11. 

За функціональним призначенням АЗОБМЗ3З поділяються на два типи: дослід- 
ницькі та спеціалізовані. Дослідницьким АЗОБМ3З характерний широкий набір функцій 
опрацювання зображень. Спеціалізовані АЗОБМ3З працюють у напівавтоматичному 
або автоматичному режимах і виконують одну визначену наперед послідовність авто- 
матичних операцій з опрацювання зображень. 

Програмне забезпечення АЗОБМ3З включає в себе такі модулі: попереднє оброб- 
лення, сегментація, вимірювання, класифікація, статистичне оброблення мікро- 
об'єктів |1-31. 

Цитологічні зображення -- це зображення окремих клітин. Мікрооб'єктами на 
цитологічних зображеннях є цитоплазма, вакуолі, ядро, ядерця, хроматин ії т.п. |41. 

При постановці діагнозу важливою проблемою є виявлення і дослідження діаг- 
ностичних ознак. Інформативною діагностичною ознакою при протіканні патологічних 
процесів є форма окремої клітини і дослідження закономірностей її зміни. На даний 
час в існуючих АЗОБМ3 відсутні модулі аналізу зміни форми клітин. Тому актуальною 
задачею є дослідження зміни форми клітин із заданою похибкою перетворення. 


Постановка задачі 


Нехай задані зображення /т; та Іт». Представимо їх у вигляді: 
р 1 2: Пр у 
Іт; «С, М/0;,, ісі2, 
де С; - контур, 0; - внутрішня область. 


Згідно з теоремою Жордана (5), для замкнутих областей існує перетворення, 
яке перетворює одну область в іншу область, тобто 0» - О(0,), де О - ідеальне пере- 


творення області О; в область 0». Реальне перетворення О" знаходиться над апрок- 
симованими областями О| і 05 за допомогою виразу 0) - О" (ОЇ). Похибку перетво- 
рення будемо визначати як площу 5 симетричної різниці областей 0, і 05, тобто 
г 
бр з 5(05АО"(О,)). 
Необхідно знайти складові похибки перетворення ду. 


Метод перетворення типу «область - область» 


Для двох заданих зображень т - С /./0, та Іт» - С» 120» області О, та О» є 
замкненими обмеженими областями на площині В). Контури Сі 1 С, є межами 
областей О, та 0», тобто С; « 20, С», « 0». Як відомо |6|, скелет області визнача- 
ється так: Др є 0; рг(Д) - р - д0| ши - р і іб|О" -р 


називають метричною проекцією точки 0. 
Для евклідової метрики використовується таке означення скелета: 


, Р" є дО. Тоді цю множину 


5К(О) з с "є дОїмножина ре(Ф) складається більше, ніж з однієї точки). 
Точка належить скелету, якщо вона є центром максимально вписаного в 
область круга. Використаємо не евклідову метрику, а /,, метрику, тобто метрику, що 


задана формулою А -8В| ш тах Їх б і, де (а, у) - координати точки А, 


бі 5 уз) - координати точки В. 
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Для /, метрики точка належить скелету, якщо вона є центром максимально 
вписаного в область квадрата зі сторонами, паралельними координатним осям. 

Як показано в роботі |71, що у випадку кусково-лінійної межі області її скелет 
також є кусково-лінійним. Крім цього доведено, що якщо область є однозв'язною, то 
її скелет є деревом. Звідси випливає, що скелет в даному випадку є скінченним об'єд- 
нанням відрізків. 

Підсумовуючи вищенаведені міркування, можна сказати, що для однозв'язних 
областей скелет є деформаційним ретрактом області, тобто 0 -» 5К(О). Скелет є 
графом без циклів, тобто його можна представити у вигляді 5К(О) - (У, Е ), де ИЙ - 
множина вершин, а Є - множина гілок. 

Перетворення типу «область - область» 0, -2 0, зведемо до перетворення 
скелетів 5К(О) ) -» 5К(О»). 

Заміна областей скелетами дає змогу перейти від двомірного до одномірного 
перетворення. Це спрощує задачу перетворення, але приводить до появи похибки 
відтворення області за її скелетом. У загальному випадку скелетизація - це некоректна 
задача |8|. Початковий скелет кожної області є зашумленим. Тому постає задача вида- 
лення шумових гілок. 

У статті |9| введено поняття є -скелету і показано, що є -скелетом називається 
скелет, для якого виконується умова: 


5К. (0) з г. є ІтО|діат рг(х)2 РА 


рг(х)з у є доЇ|х - у щ аіорно 


п 
Початковий скелет розкладається на є -скелети, тобто 5 (0) з |) 5кс, (0) Предста- 
і-1 
вимо скелет області у вигляді 5К(0) - (У, Е ,У), де - Й множина ваг гілок скелету. 
Дана множина Й - їв з Варава складається з є; ваг точок гілок е; є Е. При розкладі 
початкового скелета 5К(О) на є -скелети справедливе: 
5К(О) ЗК (0) 5 Кс» (0) 55... ЗК (0), 
тобто кожен наступний скелет є підграфом початкового графа. Ця операція відсікання 
гілок називається «обрізанням» початкового скелета. 
Нехай бк(о) - ПОХИбка відтворення області О за її скелетом. Тоді «обрізання» 
гілок скелету проводимо до виконання умови 5(5К(0), 5К, (О)) « бц(0)» Де 5 - площа 


симетричної різниці між початковою областю О та областю О", відтвореною за допо- 
могою є -скелета. З кожною точкою К є 5К(О) зв'язаний максимальний порожній 
квадрат області О(К)- (о :а(о,Ю є а(Ю)), де а(К) - сторона квадрата з центром в 
точці К. 
Об'єднання о а |) О(к) - множина максимально порожніх квадратів з цен- 
кезк(0) 
трами на гілках скелету 5К(0), яка збігається із самою областю О. 


Для оцінки похибки перетворення між областями О| -» О, використаємо симе- 
тричну різницю: дт - 5(0,40"(0,)) і бо,» бо, - похибки відтворення областей 0, 


та О» за їх скелетами. При цьому справедливе таке твердження. 
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Твердження. Верхня границя похибки перетворення між двома областями рівна: 
бт яУ бо, - бо, б 


Доведення. За означенням похибка перетворення рівна дб; - 5(0, ЛО"(0,)). 
Представимо її у вигляді бт - 0» Х0"(01) 0 0" (0| )10». Обчислимо відносну похибку 
симетричної різниці. Вона рівна: 

бт Бі бо, бог до, будо, бо, ні бо, бо, аб до, З 


Оскільки дог «« до, то бт я» бо, - бо, з 


Алгоритм визначення похибки перетворення 
типу «область -- область» 


Алгоритм визначення похибки перетворення типу «область -- область», вико- 
ристовуючи запропонований метод перетворення областей на основі перетворень 
скелетів представимо у вигляді послідовності таких кроків: 

1. Представимо контури зображень Іт; і т, в полярній системі координат 
С) -» п (0) 5 С» -» т (0). 

2. Апроксимуємо знайдені контури кусково-лінійними функціями 

п ат 
С; - )8.аувчь ве Го, ді, й дєк, і-1,2 
1 
із заданими похибками апроксимації контурів бе, , де, : 

3. Використовуючи алгоритми, які розроблені в статтях |7), |9|, знайдемо скелети 
областей Кк (О| ) 1 5К,(0») в І, метриці. 

4. Представимо отримані скелети у вигляді 50) - (Й, 5, ММ), 56(0,) - (Р, Б», 5). 

5. На основі введених множин ваг точок гілок скелетів Й та Й/, проводимо 
регуляризацію скелетів. 

6. Оцінюємо похибки відтворення скелетів областей бо, та бо, - 


7. Знаходимо похибку перетворення бу - бо, З до, - 


У свою чергу, похибка відтворення області рівна сумі похибки апроксимації 
контура області б і похибки відтворення за допомогою скелета апроксимованої 


області бно) З тобто: бо б бе а б КО) ; 


Експериментальні дослідження 


Для знаходження похибки перетворення типу «область - область» програмно 
реалізовано ряд розроблених алгоритмів: алгоритм апроксимації контуру із заданою 
похибкою апроксимації, алгоритм виділення скелету, алгоритм «обрізання» мало- 
інформативних гілок скелету на основі критеріїв втрати площі та периметру, алгоритм 
обчислення похибки перетворення типу «область -- область» та ряд інших допоміжних 
алгоритмів. Для проведення експериментів були використані цитологічні зображення, 
отримані в Тернопільському державному університеті ім. І.Я. Горбачевського. 

Для визначення функції перетворення та обчислення похибки перетворення 
типу «область - область» необхідно виконати такі кроки: завантажити вхідне цито- 
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логічне зображення та провести попереднє оброблення, за допомогою алгоритму сег- 
ментації на основі попередніх розміток виділити мікрооб'єкти на вхідному зображенні, 
визначити контурні функції мікрооб'єктів за допомогою алгоритму проходження 
контуром з можливістю зворотного ходу, провести апроксимацію контурів мікро- 
об'єктів із заданою похибкою. Приклад виділених мікрооб'єктів приведено на рис. І. 


Рисунок 1 - Приклад виділення та апроксимації контурів двох мікрооб'єктів 
ун р р тур р 


На наступному етапі здійснено обчислення координат точок, що належать ске- 
летам відповідних мікрооб'єктів. Оскільки скелетизація проводиться над мікрооб'єктами, 
представленими у цифровому форматі, отриманий скелет містить незв'язні між собою 
частини, містить «діри» та велику кількість надлишкових точок, що приводить до 
утворення великої кількості малоінформативних гілок. Приклад виділення скелетів 
наведено на рис. 2. Для усунення даного недоліку використано морфометричні опе- 
рації дилатації та ерозії, що дозволило отримати скелети мікрооб'єктів товщиною 
гілки в один піксел. 


Рисунок 2 - Приклад виділених скелетів товщиною гілки в один піксел 


Після застосування морфометричних операцій залишаються малоінформативні 
гілки скелетів, які необхідно видалити на основі критеріїв втрати площі та втрати 
периметру, що дозволило зменшити складність процесу перетворення, не перевищуючи 
при цьому допустиму похибку. 

В результаті проведення комп'ютерних експериментів було опрацьовано понад 
100 цитологічних зображень патологічних клітин долькового раку молочної залози. 

Перетворення типу «область - область» були здійснені в афінному (10) і топо- 
логічних просторах. Експерименти показали, що похибки перетворення в афінному 
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просторі в 2 - 3 рази більші за похибки перетворення в топологічному просторі. 
Максимальні похибки перетворення типу «область -- область» в топологічному просторі 
наведені в табл. 1. 


Таблиця 1 - Результати визначення похибок перетворення типу «область -- область» 


Номер Площа першого Площа другого Похибка перетворення 
експерименту мікрооб'єкта мікрооб'єкта типу «область -- 
(піксел) (піксел) область»(9/6) 
1 10832 11776 1,81 
2 8562 7456 2,03 
3 14036 12754 1,98 
4 10856 11562 2,38 
Висновки 


У роботі запропонований метод перетворення типу «область -- область» на основі 
перетворення скелетів областей. Досліджено джерела похибок при перетворенні облас- 
тей та здійснена їх оцінка. Показано, що похибка перетворення типу «область -- область» 
залежить від похибок відтворення областей зображень, які складаються з похибок 
апроксимації контурів 1 похибок відтворення областей зображень за їх скелетами. 
Проведені комп'ютерні експерименти показали, що похибки перетворення в тополо- 
гічному просторі в 2 - 3 рази менші за похибки перетворення в афінному просторі. 
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ВКЕ5СМЕ 
О.М. Вегес5ку 


Ткапз/огтапоп ої! Туре а «Кеєіоп - Кеєіоп» Еигот5 
іп те Ашотатеа Уузіетуе ої Суіїоіоєісаї Птаєе Апаїузіз 


Тре тодегп ашотагеа (0015 ої біотедіса! ітаєе ргосез5іпе, соп5і58ї плаїпіу ої 5исП 
тоаціез: подціе ої йтпаєе іприї, мубісі 15 ргезепіед а Ппебі тісгозсоре апа уїдйео сатега, 
Бу Ше рег5опа! сотригег, апа ргіогег. 50Їїмаге ої "Бе ашотагеад 10015 ої Біотедйіса! ппаєе 
ргосез8іпе, іпсіцдез пехі птодціез: ргеміоця ргосез5іпе, зертепіайоп, плеазигіпе, сіаз55ій- 
сайоп, згаці5іса! ргосез85іпо, ої ппісго обіесів. Тре ппрогапі ргобіета ої таКіпе а Фіаєпо5ія 
і5 а дДекесіїоп апа ге5еагсії ої Фабпозбіїс зуптріотя. ТРре іпбогтіпо фабпоз8іїс зупаріот аї 
Йомтіпеє ої раїбоіобіса! ргосе55е5 15 а Бог ої 5ерагаїе сеї!з апа гезеагсі ої плесрапізт ої 
її5 спапєе. Момадауз, Фе подціея ої сеїіє срапее Юогт апаїубія аге аб5епі іп Бе ехі5(їпо 
ашотаїед 10015 ої Біотедіса! ітаєе ргосе85іп2. 

Ттапзбогтаноп ої гуре «гебіоп - геєіоп» етгог5 аге соп5ідегед їп Пе агіїсіе. ТПе теїбой 
ої ітап5Їогтайоп ітаєе гебіоп5 15 ойкегей, мрбісі 15 Базей оп 58КеЇегоп8 ігап5йогтайоп ої 
тедіопя8 іп /,, плеїгіс. Тре аїсогіїит ої дегегпаїпаніоп ої ігап5богтайоп гуре «геріоп - гебіоп» 


егтог5 15 деуеіоред. Тгап5їогтаноп ої гуре «гебіоп - гебіоп» етог5 їп Бе ашіотаїей зузіетя 
ої сугоіовіса! ітаге апаїубія 15 апаЇу7ед. 

Іл Фе агіїсіе збомт, (аї ігап5їогтайоп ої гуре «геріоп - гебіоп» еггог дерепаз оп 
Бе еггог8 ої герепегаціоп ої ітаєе гебіоп5, уурісПї соп5ізі об арргохітабіоп сопіоигя еттогя 
апд еггог8 ої гесепегайоп ої ітаєе гебіоп5 оп Ба5е 5Ккеегоп5. Сотприїег ехрегітепів аге 
сопайства їп Вогіапа Дербі 7.0. епуїгоптепі. Аз а гебиії, Пе (гап5їогтацоп еггог3 п 
гороіобу зрасе їп 2 - 3 (іте5 Іе55 Бог ігап5їогтайоп еггог8 їп айпе зрасе. 


Стаття надійшла до редакції 12.04.2013. 


«Штучний інтелект» 2013 Мо 4 157 


